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We compare overlap fermions, which are chirally invariant, and Wilson twisted mass fermions in the approach 
to the chiral limit. Our quenched simulations reveal that with both formulations of lattice fermions pion masses of 
0(250 MeV) can be reached in practical simulations. Our comparison is done at a fixed lattice spacing a ~ 0.123 
fm. Several quantities are measured, such as hadron masses and pseudoscalar decay constants. 



1. INTRODUCTION 

The Wilson formulation of lattice QCD ex- 
hibits various problems: the presence of un- 
physical small eigenvalues which give rise to ex- 
ceptional configurations, the explicit breaking of 
chiral symmetry which complicates the pattern 
of operator mixing, the presence of large dis- 
cretization errors which arc reduced through the 
Symanzik improvement program. In the present 
study we consider two formulations of lattice 
QCD that are able to overcome most of these 
problems: overlap and twisted mass (tm) fer- 
mions. Overlap fermions have an exact chiral 
symmetry at finite lattice spacing a and the mass 
is an infrared cut-off, thus allowing the approach 
to the chiral limit to be performed at finite a. 
For tm fermions the twisted part of the mass also 
provides an infrared cut-off, thus solving the prac- 
tical problem of exceptional configurations which 
affects Wilson fermions. At maximal twist an- 
gle, moreover, one has automatic 0(a) improve- 
ment for various quantities like energy eigenvalues 
and matrix elements [J . The great advantages of 
overlap fermions have unfortunately the price of 
being rather expensive from the numerical point 
of view. On the other hand, tm fermions are 
rather cheap to simulate but show residual chi- 
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ral symmetry breaking effects. 

Our comparison is performed at a fixed value 
of /3 = 5.85 {a~^ ~ 1.605 GeV) and no attempt of 
a scaling analysis is performed here. Our aim 
is to investigate how both formulations behave in 
their approach to the chiral limit for a number of 
quantities like hadron masses and pseudoscalar 
decay constants. We also provide a timing esti- 
mate, from the results in ref. 

2. NUMERICAL RESULTS 

For overlap fermions we have 140 config- 
urations on 12^ X 24 lattices (L12 ^ 1.48 
fm). The bare quark masses are movO = 
0.01,0.02,0.04,0.06,0.08,0.10 and p = 1.6 0. 
The simulations for tm fermions are done at full 
twist so that all the quantities studied here ought 
to be automatically 0(a) improved [Q. In the fol- 
lowing we refer to the "twisted" basis where the 
action is the normal Wilson action with a bare 
mass term Ma- + ip-j^r^. Afcr is the bare critical 
mass determined for normal Wilson fermions and, 
for the corresponding value of the hopping pa- 
rameter, we have used = 0.16166(2) 0. The 
twisted quark mass parameter fia has been chosen 
to have the same values of the overlap masses plus 
the value 0.005 (in the plot and tables below both 
mov and /i will be called Wbaro)- We collected 140 
configurations on 12'^ x 24, 140 configurations on 
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14^ X 32 (Li4 ^ 1.72 fm) and 250 configurations 
on 16^ X 32 (Lie ^ 1-97 fm) lattices. For both, 
overlap and tm fermions, a multiple mass solver 
(MMS) has been employed. 

2.1. HADRON MASSES 

We extract hadron masses by fitting the be- 
haviour of suitable two point functions at large 
euclidean time^. We check in various ways (in- 
cluding the use of Jacobi smearing) that our 
determination is not contaminated by the pres- 
ence of excited states. Pseudoscalar masses 
are extracted from the correlation functions 
Cpi^o) = J2^{P\x)P\0)) and (only for over- 
lap) Cp_s(xo) = E^iP\x)PKO) - S''{x)S''{0)) 
where b is the flavour index and, in order to avoid 
problems of mixing with the scalar density, in the 
tm case we only consider a = 1,2. In Cp-s{xo) 
the contribution of the topological zero modes of 
the overlap operator cancels. Results are reported 
in Fig. Hand in Tab.Q 

For tm fermions we have performed simula- 
tions on three volumes and finite volume effects 
turn out to be very small for all the values 
of the mass. In the following, we will present 
only results obtained on the 16"^ x 32 lattice. 
For overlap fermions we extract the pion mass 
from Cp-s (Fig. ^ shows how large can be the 
finite volume effects due to the exact zero modes 
present in Cp, at small quark masses). The 
lowest pion mass turns out to be very small 
{M°^ ~ 230 MeV and L12 = 1-73) and thus, 
despite the cancellation above, a finite size ef- 
fect at percent level can not be excluded. In 
Fig. n] results for 0{a) improved Wilson fer- 

^In the overlap case the spectrum is automatically 0(a) 
improved. Moreover, in order to obtain 0(a) improved 
estimates of the the decay constants the bilincars can also 
be easily improved 1^. 
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Table 1 

MttO and fT^a for overlap and tm fermions. 
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Figure 1. M^a^ vs. mbarca. 

mions are also reported. In this case the simu- 
lations had to be stopped at rather large values 
of the quark mass to avoid exceptional config- 
urations. With both, tm and overlap fermions, 
we can reach instead very low values of M^. 

For overlap fermions, has, to a very good 
approximation, a linear behaviour with and 
a linear extrapolation to the chiral limit gives an 
intercept of —0.002(6). For tm fermions the be- 
haviour is better described by a quadratic form 
and the fit gives an intercept of 0.0054(4). This 
value, non-compatible with zero, is due to the 
residual 0(a) uncertainty in This uncertainty 
also induces an 0(a^/x^) effect in the pion mass. 
One can also notice from Fig. ^ that the pion 
masses obtained with tm always lay above the 
ones obtained with overlap fermions. This is due 
to the renormalization factor Z^^^ (needed to ob- 
tain the renormalization group invariant quark 
mass) that, for tm fermions, turns out to be 
roughly a factor 2.5 larger than the corresponding 
one for overlap fermions (which is close to 1 [S]). 

The vector meson mass has been extracted 
from C^(a:o) = ELi Ex (4(^)4(0)) (a-1,2) 
and C^(xo) = ELi Ex(n'(^)n'(0)) in the tm 
and overlap case respectively. The results are 
plotted in Fig.|21 In the tm case we observe (both 
with and without smearing) a progressive wors- 
ening of the plateaux for the effective masses as 
the quark mass decreases. This phenomenon is 
particularly evident for the lowest three masses 
where, due to these uncertainties, we prefer not 
to plot any result. 



S o.< 



□ 


M 

p 


overlap 


A 


M, 


twisted mass 


V 




overlap 







twisted mass 


o 




LTcupiei t^ist^"^ mass 



I 

0.1 



(Ma) 



2„2 



Figure 2. MpG and A/baryona vs. M^a 

We have also computed the proton and 
the correlators with interpolating opera- 

tors 5°==' 

^abct.,aT 



and 



e^'^^u" ^ Cjku'')ul (with fc = 1,2,3 equivalent) 
respectively. The results are presented in Fig. El 
In the overlap case, due to the smaller volume, the 
decuplet turns out to be too noisy for a reliable 
estimate. In the tm case, the same phenomenon 
as for Mp has been observed on the lowest three 
masses. Notice that, in order to obtain the phys- 
ical two point correlators, one has to rotate those 
computed in the "twisted" basis according to 
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2.2. DECAY CONSTANTS 

By using the PCAC relation, the pseudoscalar 
decay constants have been computed (without 
need of any renormalization constant) from the 
ratio f°'^ = 2mbaro|(0|P|7r)|/M2, where mbarc 
stands either for niov or for fi. Results are re- 
ported in Tab. Hand plotted in Fig. 13 At one loop 
in quenched chiral perturbation theory (qChPT) , 
/,r has neither chiral logarithms nor finite vol- 
ume effects: U / (l + (agM^/ (47r/)2) with 
/ and low energy constants. In the overlap 
case /,r nicely follows the linear behaviour pre- 
dicted by qChPT. Neglecting 317(3) breaking ef- 
fects (which are well below our statistical uncer- 
tainty) we get U = 155(11) MeV, /k = 173(8) 
MeV, /k/U = 1.11(3). 
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Figure 3. /°^a and /^'"a vs. M^a^. 

In the case of tm fcrmions, we observe a bend- 
ing of /*™ when the pion mass is small. It has 
been argued PP that for 0(a) improved quantities, 
the condition mtm ^ '^^^^qcd satisfied 
in order for the breaking of the chiral symmetry 
not to be driven by the Wilson term. The puz- 
zling fact is that, assuming a coefficient of order 
one, this condition seems to be satisfied by all of 
our data points. In Fig.O the vertical line shows 
the r.h.s of the stronger condition mtm 3> uAq^y), 
which should be valid for non-improved quanti- 
ties. We are thus left with the question of which 
inequality has to be satisfied: either the weaker 
one with a large coefficient or the stronger one 
with a coefficient of order one. This phenomenon, 
together with that observed for Mp and Mbaryon 
at the smallest quark masses, requires further in- 
vestigation at smaller values of a. 

Finally we find |2I that tm fermions are a fac- 
tor of 20-40 faster than overlap fermions and thus 
have the potential for dynamical simulations at 
realistically small quark masses on the next gen- 
eration of supercomputers. 
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